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THE EFFECT OF RADIATION ON A SMALL PARTICLE REVOLVING 

ABOUT JUPITER. 

By Theodore Henry Brown. 

1. Since the experiments of Nichols and Hull* on radiation pressure 
there has arisen the question of the effect of radiation in problems of celestial 
mechanics. Some of these problems have been already solved. Professor 
J. H. Poyntingf has worked out the conditions of size and temperature for 
two bodies in equilibrium in space and also solved to a first approximation 
the problem of two bodies moving under the pressure of radiation and the 
Newtonian attraction. Professor E. B. WilsonJ has solved the equations 
for this problem to a fourth approximation. 

The present paper treats the motion of three bodies in one plane acting 
imder radiation pressure as well as ordinary Newtonian forces: the third 
body being taken as a small satellite of Jupiter. The method of solution 
differs from the previous ones in that the time is not eliminated from the 
equations. The general effect of radiation on the particle is similar to that 
obtained in the above papers, namely that the particle is drawn in toward 
the attracting center losing eccentricity and reaching the surface after a 
finite number of revolutions. There are, however, some new effects. 
Considering the Sim to be stationary, it is clear that the particle will be in 
the shadow of Jupiter during a part of each revolution. It is foimd that 
the effect of the eclipse on the motion is to draw the particle in toward 
Jupiter and to diminish the eccentricity. But now suppose the Sim is 
moving. Then, if the motion of the particle is direct, the amount of shadow 
that the particle has to pass through will be greater than before. The 
particle will still be drawn in and lose eccentricity but the amoimt is not 
so great. On the other hand, if the motion of the particle is retrograde, 
the shadow will be less than before. Consequently the amount the particle 
is drawn in and the amount the eccentricity loses is greater. This shows 
that the number of revolutions which a particle can make with a retrograde 
motion is less than that which a similar particle can make with a direct 
motion. 

2. The assumptions in regard to the radiation pressiire will be the same 

* "The pressiire due to radiation," Astrophysical Journal, vol. 17, June, 1903. 

t Nature, vol. 75, pp. 90-93. Also "Radiation in the solar system: its effect on temperature 
and its pressiu-e on small bodies," Trans. Royal Socipty, ser. A, vol. 202 (1903). 

t "The revolution of a dark particle about a limiinous center," Annals of Math., 2d ser. vol. 
7-8 (1905-07). 
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as in Professor Wilson's paper. That is in addition to the usual gravitational 
forces there is the direct pressure of light and the Doppler reception and 
emission effects. Besides these we shall assume the motion is all in one 
plane and that the particle is started well within the sphere of action of 
Jupiter which is taken to be one twentieth the distance from Jupiter to the 
Sun. 

Let Jupiter be the origin of coordinates. Denote it by B and its mass 
by rrii. Let S be the Sun of mass mi and P the particle of mass ms. Then 
let BP = r, BS = R, SP = I, Z.PBS = <p and after Jupiter is reduced to 
rest let Bx be a line fixed in the plane of motion so that zPBx = «?, 
zSBx = a. Then, if F be the velocity of light, K the constant of radiation 
referred to S and K' that referred to B, we have for the equations of motion 
of the particle:* 






-.r^v= 



fd^Y wi2 1 K 1 dr , mi „ f r\ mi 

[Tt) = -l^-3¥vdt+T^V''^''-R)-W''''"' 

K^( r\Kldl^( r\ 



1 d f J&\ IKl d& mi„ . , K^ . 

K 1 dl„ . , mi . IK' Ird^ 

Since it is only possible to solve these equations by series and since we 
desire an approximation involving only the first power of the disturbance 
from elliptic motion, we shaE divide the equations into two parts, integrate 
each and add the results. The division is such that the first equations, 
which have a potential, represent the elliptic motion of P disturbed by the 
attraction of the Sun and the direct pressure of Ught from the Sun while the 
second set, which has no potential, gives the motion disturbed by the Doppler 
effects and the direct pressure from Jupiter. 

Take as the first set of equations 

dV fd&Y m2. rmi„/ r\ 

mi K^( r\-\ 

— -^COS (f> — yiil cos ip - ^)\ 

I d ( .d^\ [mi . mi , ^ t> • 1 

-rmV-M) = L^^'''' ^ --jT^sm ^ + jR sm ^J. 



* See E. B. Wilson, Annals of Math., 2d ser., vol. 7-8 (1905-07), pp. 135-137. 
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The right-hand sides may be expressed as partial derivatives of the function 
„ ma , TOi K mi 



(^)cos ^ + ^^i^[^(^y (3 cos^ ^ - 1) + 



mj KfrX , mi- Kri/ r' 

= -p I — 1 fiOS to -A = I - I -^r 



r 



To solve the equations we shall develop SI* in terms of the elliptic 
elements and then substitute it in the equations for the variation of the 
elliptic elements. These equations,! when the powers of e above the second 
are neglected, become 



da 2na^ dSl 






dt ~ rrii de' 






de nae dM 
dt ~ 2m2 de 


na{2 - 
2m2< 


e^)dSl 
? d-ar 



dvr ^ na{2 - e") dSl 
dt ~ 2m2e de ' 

dti _ _ 2na^ dSl nea dSt 
dt ~ mi da 2m^ de ' 

After substituting and integrating with respect to the time treating the 
elements on the right-hand sides as constants we get the first approximation 
to the changes in the elliptic elements due to the disturbing forces. These 
are, writing the first few terms: 

^ 2nc* f r3 15 ;\ cos {21 -^ 2g - 21') e , , 1 
, „ 2na' jFl , ,"|cos {V - I - g) , e .„ ., . , 



m^a mza 



say, 



* Ddaimay, M6m. de I'Acad. des Sci., vol. XXVIII (1860), pp. 33-54. 
t Tisserand, M^canique Celeste, vol. I, p. 169. 
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2mi.a 2mia 

- (m,-K) (2 - e')na' \j-cosl- .,^- ^^ „ ,. cos {Zl+2g-2l') + • 
2mi.ea'^ 1 2n 4(3n-2n') ^ ^ » -' ^^ 

" ^ 2m^a'' j - 2^ cos a - ^) - 2(^r:;^2^ cos (Z-2Z-ff) + • 
^ = ^"--^^^^W^{ 4(3ni2nO ^^^(^^ + 2, - 2^ -|^sinZ 

-4(^r^2^^^(^-2n+i«^+---} 

2m20 2TO20 

say, 

m2a'M4 ^L4 8 J 2(n - n') ^ 
2naMri , , 1 sin (r - Z - g) , 3e . ,„ , , 

+ (mx - Z) -^ {^2'} + if r^ {^5/2"}, 
2m20 2w2a 

where Z is the mean anomaly and g the longitude of perigee, the prime letters 
referring to the Sun. 

We take as the second set of equations: 

dV 









say, 

rdty dt)~ 3P V^ di Z' V dt^ ^^^ ^ 3 7 r' dt 
= T 
say. 

The solution may be obtained by expanding the right-hand sides in 
terms of the eUiptic elements and substituting in the equations for the 
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variation of elliptic elements which contain the forces on the right-hand sides. 
The expansion of dl/dt and 1/Z' in terms of r/R and cos k<p is effected as 
usual in the development of the perturbative function, the latter in turn 
being expanded in terms of the elements by use of the eUiptic developments. 
The expansions are, to the order of e desired, 



P = gi 



-o- [1 + e cos J]sin2<p — ^[1 — e cos I] sin 2tp 



— na[e sin I] cos^ <p \ — gi-^[eshil] 
+ ^^ 1 + 2e cosl + ^e^ + ^e^ cos 2l'\ - 4sr2^sm I, 
T = gi\\ — na(l + e cos I) + n'o(l — e cos Z) + -k- e sin n s ~ o ^'^^ ^*' I 

- Piyfl + c cos l\ - ff2-[l - 3e cos q, 

where ^i = KjRW, g^ = K' jZV and where v? has yet to be expanded in 
terms of the elements by use of the expression 

^ = ^- (n't-\- €'). 

Neglecting higher powers the equations to be used become 

da 2nan ^ . „ , „an 
jT = Pae sm t> + T— , 

dvr tiOi r -.-x .n ^ . \ . 1 

e -57 = — [- Pa cos t> + T{a + r) sm «>], 

(It TTli 

dei 2nfflr p 

df TO2 

In making the substitutions it will be found easier to substitute the 
expressions for P and T before expanding the functions of <p and t? into 
series of the elements. After integrating as before we get the first approxi- 
mation to the changes due to these disturbing forces. These are, writing 
only the first few terms: 

ia^ag\^-l\t+^.m2{l-V)-^.ml-lg.ot 

Am . , 2naK'e r , „„ 2t 8e . , , 
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4 L n J An' — 2n A \_n J 2(n — n) 
- sin Z — ^ sin 2e ( 1 + ^ e^ I cos Z + e cos 2Z 



2sr 
3n 






,a^.^|K-l-|cos(3Z-2Z0_3|,K_ n 
2 L «■ J 6n — 4n 2 L ^ J 



cos (Z - 2Z0 
2n - 4n' 

gri f n' ,1 J 2 grie ffie . , 

-nL«-M'°'^"3^^*-12;i"^'^-3^'"^^ 

[sin Z + e sin 2Z1 + -^ cos Z + s-V" cos 2Z, 

5.x=£cos2(Z-Z0--^cosZ--— -[z+-smZj-J^cosZ. 

The complete integral may be obtained by adding the two sets of ex- 
pressions for 5o, 5e, ehvr and 5*1 respectively. There is a constant appearing 
through the mean motion which should be added to 3o. This has been 
omitted because it does not directly concern the present problem. 

3. Consider now the effect on the motion of the particle. It is evident 
that it will suffer solar eclipse during a portion of each revolution. During 
the time it is in Jupiter's shadow the radiation pressure from the Sun will 
be cut off. We shall first determine the effect of this eclipse on the semi 
major axis and eccentricity. 

Denoting the angle which the intersection of the shadow and the orbit 
subtends at the center of Jupiter by 2/3 we have, when the particle enters 
the shadow, 

I - V = TT - ^, 

and at the time it leaves the shadow, taking the orbit to be circular, 

i _ r = IT + i3. 

Now the change in r is given by 

5r = 5o — 5e cos Z + eS-nr sin Z, 

neglecting products of the changes with e. From this we find the change 
in the value of r due to the eclipse. Since the longitude of the perijove is 
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constantly increasing, we can neglect, in making the substitution, all 
functions involving angles other than w — -bt — w'. This gives us 



Sr = Sa = r^'sin 2(Z - T) 1" ' ^sin 20. 



(1) 



Hence the particle is drawn in toward the attracting center because of the 
ecUpse. 

If for an approximation we call n a constant, we have 



P = P' 



1-^' 
n 



(2) 



where j3' is the angle between the line joining the sun and Jupiter produced, 
and the line joining the particle and Jupiter at the instant that the particle 
enters the shadow and where the particle revolves in the same direction as 
the Sim. If it revolves in the opposite direction, 



/3 = ^' 



1+^- 
n 



(3) 



If h is the time it takes the particle to make the remainder of the revolution, 
we have for direct motion 

1 



^ n n 



1 - 



and 



n 
n 



^ n 



20' 



n 



for retrograde motion. 

During this time the particle is again drawn in the additional distance 



r5 n'1 aK [5 n'l 



for direct motion and 



2t 
n 



20' 



n TT 
n 



aK 



[i-^] 



2t 
n 



20 1 

n 



(4) 



(5) 



for retrograde motion. 

In addition to this throughout the whole revolution the radiation from 
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Jupiter is affecting the particle causing it to fall in the additional distance 

AK't 



ZVan' 



(6) 



Thus the change in the semi major axis per revolution will be approxi- 
mately 



aK \ . 



2^' + 2^'^cos2^' + [f-^] 



27r - 2^' 



1-^ 

nJj 



if the motion is direct, and 
aK 



VR^ 



sm 



2S-2e'^oo.2S' + [|-f] 



2T-2/3'- 



1 + 



4XV 
ZVan 



AK'ir 
3Van 



(7) 



(8) 



if the motion is retrograde. 

Treating similarly the expression for Se we get the change in eccentricity 
per revolution. This is 



-vki^i[i^]y-^^'+^r~c^^^ 



■VL 

n 



'] 



+ 



5e n' 
4 n 



2ir - 2^' 



1 - 



_ 

n J 



7K'ex 
ZVa^n 



when the particle revolves with the sun and 



K e 

VR^l 4 



[«^rj[sin2,'-2,'^oos2,'] 



4 n . . n 



1 + 



n 



IK'eir 
SFa^n' 



when it revolves in a retrograde direction. 

Finally the particle will not make a number of revolutions iV greater 
than 



iV = 
for direct motion and 



T f 5iS:o2 + 4i!C'B2 _ ll^XoA 



ZVmi'R'' 



IT ( 5Ka^ + 4X'R^ + 5^ KaA 



(10) 
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for the retrograde motion before it reaches the siirface of Jupiter. The 
The difference of these will be 

3Vnam^ • 16- Ka^ 

N-Nj, y ~, r. (11) 

7r(5Xa2 + AK'R^) f 5Xo= + 4:K'R^ - G^XoM 

4. As a niunerical example, we shall assume the following data. 

Velocity of Ught, F = 3 X 10»«, 

Mass of the Sun, mi = 1.29 X lO^*, 

Mass of Jupiter, Wz = 1.29 X 10''', 

Distance of Jupiter from the Sun, R = 7.8 X 10". 

All of these are in the C.G.S. units.* The value of K may be determined 
from the known pressure of Ught at the earth's siurface. If the particle 
has a density about that of the earth and a radius of about .01 cm., K has 
the value 1.76 X 10"'.* We also assimie 

Constant of radiation of Jupiter, K' = K X 10~*, 

Distance of the particle from Jupiter at an apse, a = ^R X 10"''. 

Then if we suppose the shadow cast by Jupiter is a cylinder whose base 
is a great circle of Jupiter, j8' will be equal to .0178 radians or a Uttle more 
than 1°. Then it will be found from (1) with (2) or (3) that, due to the 
eclipse, the direct particle will fall in a distance 8.98 cm. while the retrograde 
particle will fall in a distance 9.40 cm. It may be Jioted that these values 
increase to a maximum when the particle has fallen in so that ^' equals 
ir/4 and then decreases to nearly when the particle skims the surface of 
Jupiter. During the rest of the revolution equations (4), (5) and (6) show 
that the direct particle falls in a distance .688 km., while the retrograde 
particle falls ta a distance .689 km. These will be approximately the totals 
as the ecUpse effect is not quite appreciable at this distance. 

Treating the expression for the change in the eccentricity in the same 
way we find the ecUpse change for the direct particle is .122 X 10~^''e while 
that for the retrograde particle is .117 X 10~^''e. During the rest of the 
revolution the change for the direct particle will be 19.15 X lO^e which is 
approximately that for the retrograde particle. Thus the total change will 
be 19.16 X 10-»e for both particles. 

To find the greatest number of revolutions we will start the particle a 



* Wilson, Annals of Math., 2d ser., vol. 7-8 (1905-07), p. 148. 
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little farther out from Jupiter and let 

Distance from Jupiter, a — R X 10~^. 

Then the niunber of revolutions which is found from (9) or (10) probably 
will not be more than 5.87 X 10' for the direct and 5.54 X 10' for the 
retrograde particle. The difference between these, given by (1 1) , is 3.3 X 10*. 
This indicates that the retrograde particle is drawn in much faster, a fact 
which can be verified in general by comparing equations (7) and (8). 



